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The Unitary Fermi Gas (UFG) is one of the most strongly interacting systems known to date, as 
it saturates the unitarity bound on the quantum mechanical scattering cross section. The UFG 
corresponds to a two-component Fermi gas in the limit of short interaction range and large scat- 
tering length, and is currently realized in ultracold-atom experiments via Feshbach resonances. 
While easy to define, the UFG poses a challenging quantum many-body problem, as it lacks any 
characteristic scale other than the density. As a consequence, accurate quantitative predictions of 
the thermodynamic properties of the UFG require Monte Carlo calculations. However, significant 
progress has also been made with purely analytical methods. Notably, in 2005 Tan derived a set 
of exact thermodynamic relations in which a universal quantity known as the "contact" C plays 
a crucial role. Recently, C has also been found to determine the prefactor of the high- frequency 
power-law decay of correlators as well as the right-hand-sides of shear- and bulk viscosity sum 
rules. The contact is therefore a central piece of information on the UFG in equilibrium as well 
as away from equilibrium. In this talk we describe some of the known aspects of Fermi gases 
at and around unitarity, show our latest Monte Carlo results for the contact at finite temperature, 
and summarize the open questions in the field, some of which we are starting to answer using 
large-scale Monte Carlo calculations by adapting methods from Lattice QCD. 
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1. Introduction 



In recent years, our understanding of universality [jl|] in non-relativistic many-body quantum 
mechanics has increased dramatically. By universality we mean independence from the details 
of the interaction, in the same sense as in the context of second-order phase transitions, but in the 
absence of long-range correlations throughout the system (except potentially at special points in the 
phase diagram). Among the systems displaying this property, perhaps the most dramatic example 
is the so-called unitary Fermi gas (UFG). This system is a two-component Fermi gas tuned to the 
limit of vanishing interaction range r and large s-wave scattering length a s , and is termed "unitary" 
because it saturates the unitarity bound imposed on the scattering cross section by the unitarity of 
quantum mechanics. In short, the unitary gas is a resonant quantum mechanical many-body system. 

A few years ago, the unitary limit was realized in metastable ultracold atomic clouds in various 
laboratories around the world [Q] and it has been under intense scrutiny by the atomic, molecular 
and optical physics (AMO) community ever since |^|j. Interest in the UFG transcends those ar- 
eas, however, with a considerable amount of research being carried out within the nuclear physics 
community well before and after the first AMO experiments [Q]. This, of course, is itself a mani- 
festation of the universality of the UFG, as nuclear systems characteristically display short ranges 
and unnaturally large scattering lengths, although in a vastly different absolute scale than atomic 
clouds, the natural scale being in each case the Fermi momentum k F . 

More recently, the limit t f r < 1< k F a s — > °° has been shown to imply non-relativistic 
conformal invariance, as described in Ref. [S]. In turn, this results in a set of non-trivial relations 
between the system in homogeneous space and in a harmonic trap, as first shown in Ref. [|6L[7p. In 
a separate line of research, short-distance correlations were shown by Tan [||] and others [Q, [T(J to 
be completely encoded in a quantity C, which Tan called the "contact". Specifically, we may define 
the contact as 

C=]imk 4 n a (k), (1.1) 

where n a (k) is the momentum distribution for spin a expressed as a thermal average. This re- 
markable property stems in part from the short-range nature of the interaction, which implies that 
at resonance the many-body wavefunction is essentially that of a free gas, with the added bound- 
ary condition that it diverges as 1/r when two coordinates are set to a short distance apart r JTT|]. 
Following the work of Tan and others, the last couple of years have seen considerable activity 
extending the analysis of short-range correlations in many-body systems to systems away from 



unitarity as well as to different dimensions [ |12| ] and to a growing set of thermodynamic and even 



hydrodynamic quantities Q13| , |14Q . The latter, in particular, point to the fact that C is relevant not 



only in equilibrium but also away from equilibrium (see Ref. [ |15| ] for a comprehensive review). 

In spite of experimental advances and progress from the formal and analytic points of view, 
the UFG remains a challenging many -body problem. The reason for this is that, while resonant and 
therefore strongly interacting, the UFG has as few scales as a non-interacting Fermi gas. Enhanced 
symmetries aside, such a lack of scales implies lack of small parameters to perform an expansion, 
such that non-perturbative numerical methods are required. Indeed, while we know that Tan's 
contact plays a crucial role in the dynamics of the UFG, the only way to determine it accurately 
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and reliably is by using numerical methods such as Quantum Monte Carlo (in any of its various 
incarnations, in particular on the lattice). 

The work shown here represents the first attempt to determine the contact at finite temperature 
in a non-perturbative fashion. To this end, we have adapted methods from Lattice QCD, namely 



Hybrid Monte Carlo [16] and applied them to the calculation of the momentum distribution of the 



UFG. The next section outlines the main features of the algorithm and Sec. || shows our results and 



conclusions, which were first published in Ref. [18]. 



2. Algorithm & lattice formulation 



The lattice formulation we have used for this work follows closely that of Ref. [|19|], but differs 
in at least three notable aspects. Firstly, we determine the bare lattice coupling constant g corre- 
sponding to the unitary regime by using Luscher's formula pOO as in Ref. [||], without imposing 



a spherically symmetric cutoff. This procedure yields g ~ 5.144 in the unitary limit. Secondly, we 



use the compact, continuous Hubbard-Stratonovich J22| ] transformation 

1 t K 

exp (Tgn t ,-%) = — J_da i [l +Bsin(a,-)%] [l +Bsin(a,-)%] , (2.1) 

where a, (not to be confused with the spin projection) is the auxiliary field, with B 2 /2 = exp(rg) — 
1, and X denotes the lattice spacing in the imaginary time direction. We find that a time step x ~ 0.05 
is sufficiently small to render temporal discretization errors insignificant. The above representation 



(referred to as "Type 4" in Ref. [23]) was found to be superior with respect to acceptance rate, 



decorrelation and signal-to-noise properties than the more conventional unbounded and discrete 



forms [24]. Finally, the use of a continuous auxiliary field allows us to perform global updates 



using the Hybrid Monte Carlo (HMC) algorithm Q16|]. Our implementation of the HMC algorithm 



does not use pseudofermions but rather relies on a direct calculation of the fermion determinant in 



a purely spatial rather than spacetime formulation Q17|]. In addition, we use Fourier acceleration 
to propagate states in imaginary time. This enables global updates at all temperatures and lattice 
sizes, and scales approximately as ~ V 2 logV (at fixed temperature) for moderate spatial lattice 
volumes V, to be contrasted with the ~ V 3 scaling of approaches based on local updates. 

3. Results and conclusions 

We have performed calculations at zero as well as finite temperature, in the former case using 



an approach similar to Ref. [|23]]. Our main results correspond to 40 — 50 particles at N x = 10 
and 70 — 80 particles at N x = 12, in addition to limited data for N x = 14. In Fig. |l| (left panel), 
we show the momentum distribution n(k) as a function of temperature T/s F . We have computed 
n(k) by averaging over the angular directions on the lattice as well as over the imaginary-time 
slices. In this way, we find that ~ 200 uncorrected auxiliary field samples for each datapoint gives 
excellent statistics for n(k). Multiplying n(k) by k 4 , as plotted in the right panel of Fig. [IJ we find 
a maximum at k ~ k F and a leveling out at high momenta, with the asymptotic regime setting in 
at approximately 2k F , at the lowest temperatures. There is no a priori reason for the asymptotic 
regime to set in at such low momenta; our work is the first to point out this fortunate situation. 
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We study the temperature dependence of this "plateau", which allows us to determine the contact 
C/ {Nk F ) as a function of temperature. The corresponding results are given in Fig. ||[ together with 
a comparison with other theoretical predictions. Our results indicate that n(k) follows the expected 
~ k~ 4 dependence accurately up to at least k ~ 4k F , at which point the signal deteriorates, possibly 
due to lattice effects. 
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Figure 1: (Color online) Right panel: Momentum distribution n(k) from QMC for N x = 10 as a function 
of k/k F , for various temperatures ranging from zero to T/e F ~ 0.5. The solid lines are intended to guide 
the eye, and the statistical errors are of the size of the symbols. Inset: n(k) for N x = 14 in a log-log scale, 
showing the asymptotic ~ k~ 4 behavior. Left panel: Plot of 3n 2 (k/k F ) 4 n(k) for N x = 12 as a function of 
k/kp at T /e F = 0.178 and 0.404. The "plateaux" at large k/k F give the intensive dimensionless quantity 
C/(Nk F ). At low T/e F , the asymptotic region is reached at k/k F ~ 2. Inset: = 10 results at T = 
showing only slight dependence on the time step t. 



Our results show that the contact C grows with temperature well beyond the superfluid phase, 
which is suggestive of a peak C max ~ 3.4 at T/e F ~ 0.4. This scenario agrees qualitatively with 
Ref. [25], as well as Ref. [|26|]. Since C measures the number of particle pairs (of both spins) whose 
separation is small, the appearance of a maximum indicates an enhancement in such short-range 
correlations. We find the scale at which the k~ 4 law sets in (see Fig. [I|) to be k ~ 2k F at finite T/e F 
and somewhat lower for the ground state. 

In summary, we have computed the momentum distribution n(k) and the contact Cj (Nk F ) for 
the UFG at zero and finite temperature, using a lattice formulation of the many-body problem, in 
conjunction with the HMC algorithm. Our results represent the first fully non-perturbative calcula- 
tion of n{k) free of uncontrolled approximations. We find that the contact at T = takes the value 
~ 2.95 ±0.10 and increases as a function ofT/e F in the low- and intermediate-temperature regimes 
that we have explored, which is consistent with the phonon-dominated scenario of Ref. [25]. Our 
results complement the calculations of Refs. [25, 28], and are suggestive of a maximum 



in C/(Nk F ) at T/e F ~ 0.4, which agrees qualitatively with Ref. [ |26| ] but disagrees with Ref. [[27|]. 
While calculations at higher temperature T /e F ~ 1 are feasible, an improved understanding of the 
finite density effects is clearly called for. 
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Figure 2: (Color online) Summary of QMC results for C/(Nk F ) as a function of T/e F , as determined 
from the large kjk F behavior of n{k). The errorbars are dominated by systematics related to the residual 
fluctuations in the plateaux, as shown in the previous figure (right panel). Also shown are the t-matrix 
calculations of Ref. [ p6[ , |2^ ], the virial expansion of Ref. [^8| and the diagrammatic Monte Carlo result of 
Ref. [|§]. 
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